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Abstract The mirror dark matter (MDM) model of 
Berezhiani et al. has been shown to reproduce observed 
galactic rotational curves for a variety of spiral galax- 
ies, and has been presented as an alternative to cold 
dark matter (CDM) models. We investigate possible 
additional tests involving the properties of stellar or- 
bits, which may be used to discriminate between the 
two models. We demonstrate that in MDM and CDM 
models fitted equally well to a galactic rotational curve, 
one generally expects predictable differences in escape 
speeds from the disc. 

The recent radial velocity (RAVE) survey of the 
Milky Way has pinned down the escape speed from the 
solar neighbourhood to v csc = 544^4g km s _1 , placing 
an additional constraint on dark matter models. We 
have constructed an MDM model for the Milky Way 
based on its rotational curve, and find an escape speed 
that is just consistent with the observed value given 
the current errors, which lends credence to the viabil- 
ity of the MDM model. The Gaia-ESO spectroscopic 
survey is expected to lead to an even more precise es- 
timate of the escape speed that will further constrain 
dark matter models. However, the largest differences 
in stellar escape speeds between both models are pre- 
dicted for dark matter dominated dwarf galaxies such 
as DDO 154, and kinematical studies of such galaxies 
could prove key in establishing, or abolishing, the va- 
lidity of the MDM model. 

1 Introduction 

Over the last decade, mirror matter has emerged as 
a promising candidate for dark matter. The idea that 
there might be a hidden "mirror sector" of the Universe 
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dates back to [2JJ , was subsequently expanded upon by 
[22] , and owes its modern form (in the context of gauge 
theories) to [18]. It is based on the observation that 
parity symmetry is violated by the weak interaction in 
the "ordinary" sector we observe, but may be restored 
by the existence of a duplicate sector with the same 
fundamental particles and microphysics (apart from the 
opposite handedness of the weak interaction) . Support 
for the existence of a mirror world may also come from 
the recent experiments of [B] on neutron oscillations. 

Mirror dark matter (MDM) models linked to su- 
perstring inspired unifications of the four forces in the 
early universe have been discussed by [23) . The ordi- 
nary and mirror sectors are coupled through gravity 
but not the other fundamental interactions, which ex- 
plains the intrinsic suitability of mirror matter as a 
dark matter candidate. Much research has been done 
on cosmological signatures of mirror matter; these turn 
up in the context of early Universe thermodynamics, 
Big Bang nucleosynthesis, primordial structure forma- 
tion/evolution, cosmic microwave background (CMB) 
radiation and large-scale structure power spectra. De- 
tailed reviews of this work have been provided by [3] 
and QI]. 

Astrophysical signatures of mirror matter have also 
been looked at in the context of galactic rotational 
curves, using a specific model for mirror gravity de- 
veloped by [32], @] and [5] (henceforth BPR). In this 
model, the ordinary and mirror sectors are assumed to 
have separate metric tensors (which are coupled to a 
third metric tensor), leading to a Yukawa- like modi- 
fication of the gravitational interaction between both 
sectors. An ordinary matter test particle located at a 
radial distance r from a point-like source of ordinary 
mass Mi and mirror mass M2 will feel a Yukawa-like 
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potential 

(r) = ~ (Mi + M 2 + (Mi - M a ) e~^) , (1) 

where G is the Newtonian gravitational constant and 
r m is the Yukawa radius. The free-fall acceleration of 
the particle follows as 



9(r) 



G 
2^ 



Mi + M 2 + (Mi - M a ) 1 



(2) 



Although mirror matter is unable to form extended 
halos around galaxies due to its collisional and dissipa- 
tive nature, the combined model of mirror matter and 
mirror gravity has been successfully fitted to the ob- 
served rotational curves of several disc galaxies by BPR, 
on the assumption that ordinary matter has a similar 
density profile to its mirror counterpart. This makes it 
the most viable theory of MDM that is presently avail- 
able. 

In this paper, we investigate some implications of 
the BPR model over and above its demonstrated suc- 
cess in modelling rotation curves, and suggest astro- 
physical signatures that may be observable and poten- 
tially used to discriminate between the MDM and CDM 
models. In Section 2, we use point-like source models 
to investigate the orbits of ordinary matter particles in 
the presence of mirror matter, and show that significant 
differences are expected in the nature of orbits near the 
Yukawa radius. We then construct in Section 3 CDM 
models for the galaxies NGC 2403, DDO 170 and DDO 
154 (whose rotational curves have previously been mod- 
elled using MDM), and calculate and compare escape 
speeds from the stellar disc as a function of galactic ra- 
dius for both types of model. We show that the MDM 
models predict specific values for the escape speeds that 
are generally different from those in the CDM models. 
In Section 4, we construct MDM and CDM models for 
the Milky Way, and find that the escape speed from 
the solar neighbourhood in the MDM model is in good 
agreement with the current observed value. We discuss 
and present our main conclusions in Section 5. 



2 Point-like source orbits 

The deviation of mirror gravity from an inverse-square 
law leads to distinctive orbital dynamics in the gravita- 
tional field of a point-like source, especially at distances 
r ~ T m - Two prominent features are orbital precession 
and a varying escape speed to circular speed ratio, nei- 
ther of which are present for point-like sources in clas- 
sical gravity. 



Orbits for an ordinary test body in the gravitational 
field given by Equations (HJ and @ are described by 
its Binet equation 
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Mi + M 2 + (Mi - M 2 ) 1 + 



(3) 



where the reciprocal radial distance u = 1/r is a func- 
tion of the plane polar angle 8, and the specific angular 



momentum h = L/r 



'9 is a conserved quantity. 



Taking the Yukawa radius r m — 10 kpc (as assumed 
by BPR) here and henceforth, we solve Equation (3) 
numerically for orbits about a point-like source with 
Ma/Mi = 10 (see Figure H}. 

The test body is given an initial radial distance of 
ro = ar m and an initial transverse velocity of vo = 
J GMi/tqO, which is the circular velocity in the ab- 
sence of mirror matter and mirror gravity. Its orbit is 
essentially Newtonian for a -C 1, as gravitational inter- 
action between ordinary and mirror matter is negligible 
at small distances, and Equation @ reduces to 



GAh 



in the limit r — > 0. 

For a 3> 1, gravitational interaction between ordi- 
nary and mirror matter becomes universal, and Equa- 
tion ([2]) reduces to 



9(r) = 



G (Mi + Ma) 
2r 2 



in the limit r — > oo. The test body feels the combined 
inverse-square field of both matter types (but with an 
effective gravitational constant of G/2), resulting in an 
elliptical orbit. 

Precessing orbits are manifest for a ~ 0.1 to a ~ 10, 
and are explicitly non-classical. They may be expected 
to play a role in galactic stellar dynamics at distance 
scales on the order of the Yukawa radius from the centre 
of a galaxy. 

A more useful signature of mirror gravity, which 
may allow it to be observationally corroborated or fal- 
sified, lies in the escape speeds that it predicts. For a 
point- like source with gravitational field given by Equa- 
tions (H|) and ([2]), the escape and circular speeds are 
given respectively by 



v csc (r) = yj—2<f> (r) 



= \/~ (M, + M> - (Mj - M 2 )< 



(4) 




Fig. 1 Mirror gravity orbits (solid) for ordinary test body about point- like source with M2 /Mi = 10, initial radius ro = ar m 
and initial velocity Vq = ^GMi/rgO. The corresponding Newtonian orbits (dashed) are circular. All orbits are shown up to 
five revolutions. 
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Fig. 2 Plot of v eBC /v C it against r/r m for M 2 /Mi = 1000 
(solid), M 2 /Mi = 100 (dashed) and M 2 /M 1 = 10 (dot- 
dashed). The Newtonian ratio i> eB cA>cir = y/2 is recovered 
for M 2 /Mi = 1 (dotted). 



and 



■ (r) = \J rg (r) 



Mi + M 2 + {Mi - M a ) 1 + 



(5) 



In general, the ratio v esc /v c i x ^ \[2 and is a func- 
tion of r, in contrast to the Newtonian case. A plot of 
Wcsc/^cir for different values of M2/M1 > 1 shows that 
Wcsc/^cir > v2) with the discrepancy being more pro- 
nounced for larger mass ratios (see Figure [2]). Q This 
suggests that mirror gravity may predict significantly 
larger escape speeds compared to classical gravity when 
both models are fitted to give similar circular speeds, 
as is done in the modelling of galactic rotational curves. 



3 Escape speeds in disc galaxies 

Mirror matter and mirror gravity have been put forth 
by BPR as a collisional and dissipative dark matter 
candidate capable of explaining flat galactic rotational 
curves. They have successfully fitted this MDM model 
to the observed rotational curves of several disc galax- 
ies, using simplified and similar mass density profiles 
for the visible and dark matter distributions. In this 
section, we obtain comparable fits using a CDM model 
of the same complexity, and generate the corresponding 
galactic escape speed curves for both models. 



1 The cosmological ratio between the dark and visible energy 
densities is P2/P1 ~ 10 in the BPR model. M2/M1 should ap- 
proach the cosmological ratio for larger galaxies, while smaller 
galaxies are expected to be more dark matter-dominated. In 
any case, only values of M2/M1 > 1 will be considered in this 
paper. 



The distribution of visible matter in a typical disc 
galaxy is approximated most simply as a two-dimensional 
disc with the exponential surface mass density 



A,d (r) 



M hd 

?7rr- 2 



(6) 



where M\^ and ti^ are the total mass and characteris- 
tic length scale of the disc respectively. In the BPR 
model, the distribution of mirror dark matter is as- 
sumed to take on a similar profile (i.e. another expo- 
nential disc concentric to the visible disc); this is jus- 
tifiable since mirror matter has the same microphysics 
as ordinary matter and can be expected to undergo the 
same process of galaxy formation, up to a rescaling of 
parameters. Hence its surface mass density is given by 



Z2,d (r) 



2.d 



(7) 



where Mi^ and r-i^d are the corresponding parameters 
for the mirror dark disc. 

For an ordinary test body with position vector r 
in the plane of the two discs, the acceleration vector 
follows as 



g 0) = - 



G 



d 2 v' 



r — r 



+ (S hd (r')-S 2 , d (/)) ! + 



2i,d (/) + S 2M (V) 
|r-r'|\ 



(8) 

where g (r) = — g (r) f . The integration can be per- 
formed semi- analytically in this case, and yields the 
rotational curve v (r) = y/ rg (?*) of the model. 

Assuming that the mass-to-light ratio 7 = Mi^/Ml 
of the visible matter (where the luminous mass Ml 
is an observed fixed parameter) is constant for each 
galaxy, the dimensionless quantities 7, ft = Mi^jMi^ 
and a = ri^jfiA have been taken by BPR as fit pa- 
rameters for the rotational curves of four disc galaxies. d 
The visible length scale 7*1 <j for each galaxy is deduced 
from and fixed by observations. The resulting three- 
parameter fits are demonstrably good, with Xrcd ~ 1 
for each galaxy (see Figure [3] and Table [T]) . 

Numerous galactic rotational curve fits for CDM 
models exist and are well documented in the litera- 
ture. These are of varying complexity, but nearly all of 
them incorporate components such as a galactic bulge 
or a gaseous contribution, in addition to the standard 



2 Rotational curve fits for the large spiral galaxy NGC 7331 
require the contribution of an additional bulge component, 
and will be omitted from the MDM-CDM(IT) comparisons 
in this section. We address the issue of a galactic bulge when 
modelling the Milky Way in Section U 
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galactic disc and dark halo components. Indeed, the ro- 
tational curves of the galaxies considered by BPR have 
been fitted for a multi-component CDM model by [2] 
(henceforth BBS). 

We adopt a simplified two-component version of the 
BBS model in this paper, allowing direct comparison 
between the MDM model and a CDM model of the 
same complexity. The distribution of visible matter is 
again an exponential disc with surface mass density 
given by Equation ©, while a pseudo- isothermal trun- 
cated sphere (which we shall refer to from now on as 
CDM (IT)) is assumed for the distribution of dark mat- 
ter. 

The volume mass density of the dark halo is given 

by 



Ph (r) 



1 + 5 



(9) 



Furthermore, we have generated synthetic rotational 
curves for two disc galaxies, one whose physics is dic- 
tated by MDM and the other by CDM(IT). Typical ob- 
servational errors (Gaussian noise) have been assigned 
to these synthetic galaxies, which are then treated as 
actual observed galaxies and fitted using the procedures 
outlined in this section. We find that the MDM fit re- 
covers best-fit values in good agreement with the ac- 
tual values for the synthetic MDM galaxy, while the 
CDM(IT) fit matches the parameters of the synthetic 
CDM(IT) galaxy as well; this verifies the accuracy of 
our fitting algorithms. Finally, a side-by-side compari- 
son of all the MDM and CDM (IT) fits shows that both 
models are able to reproduce both the observed and 
synthetic galactic rotational curves, and equally well 
(see Figure [3] and Tabled]). 



where p s and r s are the central density and character- 
istic length scale of the halo respectively. This profile 
avoids the central cusps inherent in those produced by 
numerical n-body simulations (e.g. the model by |31jV 
but has a divergent total mass and must be truncated 
at a suitable distance from the galactic centre. 

An analytical expression is available for the rota- 
tional curve of this simplified BBS model. We can write 



V M = \Adisc W+"halo 



(10) 



where i>di sc ( r ) an d Vhnio (r) are rotational curves for the 
visible disc and dark halo respectively. The contribution 
from the disc has been obtained by [TO] , using a method 
by [35], as 



w disc 



(r) 



GM hd r 2 



2r M 



Io 



2r M 
r 



K Q 



2ri,d 



2r M 



(11) 



where the /„ and K n are modified Bessel functions eval- 
uated at r/2ri :( i. [I] gives the contribution from the halo 

as 



J halo 



(r) 



1 - — tan" 1 — 



(12) 



where w max = A^Gp s r\ is the asymptotic rotational 
speed for the halo at r 3> r s . 

As in the BPR model, we use data from [T] , [TJ] , [25] 
and [5] for the intermediate spiral galaxy NGC 2403 
and the dwarf galaxies DDO 170 and DDO 154. Again 
assuming a constant mass-to-light ratio for each galaxy 
and taking Ml and ri lC j as fixed parameters, we fit the 
rotational curves of these three disc galaxies by using 
7, Vmux and r s as fit parameters. 



Now, for the MDM model, the gravitational poten- 
tial felt by an ordinary test body with position vector 
r in the plane of the two discs is given by 



0MDM (r) 



1 



disc l r — r I 

'I" 



+ (27 M (/) - S 2A (r')) e 



(13) 



where i^mdm (?) depends only on r due to cylindrical 
symmetry. For the CDM(IT) model, the gravitational 
potential felt by a test body with position vector r in 
the plane of the disc is given by 



^cdm(it) (r) = -G 



<Fv> ElAr ') 



-G 



disc 



halo 



(14) 



where <^cdm(it) ( r ) depends only on r as well. 

Using the parameters given in Table [T] we obtain 
MDM and CDM (IT) escape speed curves v (r) = ^-2<j) (r) 
for the five disc galaxies by numerical integration, trun- 
cating the discs and halos at a value of r just beyond 
the farthest data point for each galaxy (see Figure 2]). 

As anticipated from our discussion of point-like source 
models in Section^ escape speeds for stars in the galac- 
tic plane are mostly higher in the MDM model. The 
offset between the two models appears to be nearly 
constant out to the truncation distance. It is both abso- 
lutely and relatively larger for the dwarf galaxies (Av > 
30 kms" 1 and Av/v CDM{1T) > 0.2) than for NGC 2403 
(Av ps 10 kms -1 and Z\v/«cdm(it) ~ 0.03). 

This systematic escape speed discrepancy may be 
attributed to the increasing importance of mirror dark 
matter over ordinary matter in the smaller galaxies (as 



(i 



NGC 2403 



DDO 170 




Fig. 3 Observed and synthetic rotational curves for different disc galaxies, with both MDM (solid) and CDM(IT) (dashed) 
three-parameter fits. Parameters and values of Xred are gi ven m Table [TJ 



indicated by the higher values of which in turn is 
due to typical distance scales in these galaxies being 
smaller than the Yukawa scale r m (~ lOkpc) such that 
dark matter is less effective in gravitationally binding 
to ordinary matter. In the largest galaxies such as the 
Milky Way (see Section 4), the ratio of the energy den- 
sity of mirror dark matter to visible matter P2I P\ ap- 



proaches the cosmological value 
BPR model. 



10 as required by the 



Crucially, we have also observed in our analysis that 
escape speeds for the MDM model do not depend signif- 
icantly on the choice of truncation radius, and are con- 
strained by the rotational curve fits. Thus, observations 
of the stellar escape speeds in disc galaxies of different 
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Table 1 MDM and CDM(IT) best-fit parameters for the rotational curves of different disc galaxies. Ml and r\ jd are taken 
as fixed parameters, while 7, /?, a, w max and r a are taken as fit parameters. Both models provide realistic (all values of 7 ~ 1) 
and good (all values of Xrcd ~ -0 n ^ s - ^ ne fitting algorithms recover best-fit values in good agreement with the actual values 
for the synthetic disc galaxies (given in parentheses). 







NGC 2403 


DDO 170 


DDO 154 


Synthetic (MDM) 


Synthetic (CDM(IT)) 


Fixed parameters 


M L (10 9 M G ) 


7.9 


0.18 


0.05 


8.0 


0.10 




n,d (kpc) 


2.0 


1.1 


0.5 


1.6 


0.8 


MDM fit 


7 = M ltd /M L 


1.6 


2.2 


0.6 


1.0 (1.0) 


1.2 




M 1)d (lO 9 M ) 


12.6 


0.4 


0.03 


8.2 (8.0) 


0.12 




j3 = M 2id /M 1<d 


20.6 


190 


2000 


22.3 (23.0) 


567 




M 2 , d (lO 9 M ) 


260 


76 


60 


183 (184) 


68 




a = r 2 ,d/ri,d 


1.6 


0.8 


1.3 


1.6 (1.6) 


1.0 




r2,d (kpc) 


3.2 


0.87 


0.66 


2.5 (2.6) 


0.8 




•^red 


1.3 


1.5 


0.6 


1.3 


1.3 


CDM(IT) fit 


7 = M 1 . d /M L 


1.5 


1.1 


0.2 


1.0 


0.9 (1.0) 




M!, d (10 9 M o ) 


12.1 


0.2 


0.01 


7.7 


0.09 (0.10) 




»mK (kms -1 ) 


155 


77 


61 


129 


70 (70) 




r s (kpc) 


5.0 


2.1 


1.7 


4.1 


2.0 (2.0) 




V 2 
^rcd 


1.6 


1.2 


0.8 


1.3 


1.0 



sizes could prove key in establishing — or abolishing- 
the validity of the MDM model. 



We employ the tractable Plummer sphere for this 
purpose, in lieu of the truncated power-law bulges used 
in more detailed mass models. Hence the gravitational 
potential of the visible bulge is given by 



4 Escape speeds in the Milky Way 

Although the rotational curves of numerous external 
disc galaxies have been obtained to a good degree of 
precision, there are no direct estimates of their escape 
speeds at the present time. Conversely, while the rota- 
tional curve of our own galaxy is less established, ef- 
forts have been made by pH] to constrain the escape 
speed at the solar neighbourhood through the Radial 
Velocity Experiment (RAVE) survey. The current esti- 
mate is significantly in excess of v^^cir (where i> c i r ~ 
220 kms -1 is the local circular speed), and has been 
taken as independent evidence in favour of CDM mod- 
els. In this section, we repeat the MDM-CDM(IT) com- 
parison for the Milky Way, and show that the MDM 
model is also able to predict a local escape speed com- 
parable to the currently observed value. 

As in Section[31 the visible disc of the Galaxy in both 
models is taken as an exponential disc and described by 
Equation ©. Likewise, the density profiles of the mirror 
dark disc in the MDM model and the dark halo in the 
CDM (IT) model are given respectively by Equations 
([7]) and In both models, however, an additional 
bulge component is required to reproduce the observed 
Galactic rotational curve. 



'/'bulge (r) 



GM 1 



(15) 



1,6 



where Mi & and r^j, are the total mass and character- 
istic length scale of the bulge respectively. Its volume 
mass density follows as 



Pi,b (r) 



3A/i, 6 
47rrf . 



r 



(16) 



from Poisson's equation. 

By assumption, the presence of a visible bulge neces- 
sitates the addition of a mirror dark bulge component 
to the MDM model. This is similarly described by 



P2,b 0) = 



3M- 



2, ft 



r 



(17) 



where M^.b and r^.b are the corresponding parameters 
for the dark bulge. The acceleration vector of an ordi- 
nary test body with position vector r in the Galactic 
plane follows as 



g 0) = gdisc (r) + gbuige (r) 



(18) 
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Fig. 4 Predicted MDM (solid) and CDM(IT) (dashed) escape speed curves for stars in the galactic plane of observed and 
synthetic disc galaxies, using parameters given in[T] The actual (dotted) escape speed curves for the synthetic galaxies have 
been included as well. 



where gdisc (r) is given by Equation flSJ) and gbuige ( r ) is 
obtained in identical fashion by numerically integrating 
Equations (flo]) and (fT7|) over the bulge. 

With the addition of the visible bulge, the rotational 
curve of the CDM(IT) model becomes 

« M = V^disc ( r )+<io W + ^uigo W. (19) 

where v^ isc (r) and w^aio ( r ) are gi ven by Equations 
and dm as before, and v^ ul (r) follows from Equation 



dD as 

2 d GAh.y ( r 2 V § 

<W (r) = r—fani gc (r) = ^ I rf7 ) ' 

(20) 

For the MDM model, we assume a constant dark- 
to-visible mass ratio (such that (3 — M2,i/M-yj) and 
a constant dark-to-visible length scale ratio (such that 
a = T2.il r\.i) in the Galaxy. We use M\^i ^i,bi P and 
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Fig. 5 Observed Galactic rotational curve out to 2Rq , with 
both MDM (solid) and CDM(IT) (dashed) four-parameter 
fits. Parameters and values of x? e d are gi ven m Table 2. 



Table 2 MDM and CDM(IT) best-fit parameters for the 
Galactic rotational curve out to 2Rq. T\ r cL an <A ?%(, are taken 
as fixed parameters, while Mi^, Mi j, ft, a, w max and r s are 



taken as fit parameters. Both models provide good (Xred 
fits that agree with observation and other mass models. 



1) 





Parameters 


Values 


Fixed parameters 


r%,d (kpc) 


3 




ri,6 (kpc) 


0.3 


MDM fit 


Af lid (lO 9 M ) 


48 




M ub (10 9 M Q ) 


13 




P = M 2 ,i/M!,i 


12 




M 2]d (10 9 M Q ) 


576 




M 2 , b (10 9 M Q ) 


156 




a = r 2 ,i/ri :i 


1.7 




r 2 ,d (kpc) 


5.1 




r 2 ,b (kpc) 


0.51 




<^red 


1.3 


CDM(IT) fit 


M 1 , d (10 9 M Q ) 


54 




Mi, b (1O 9 M ) 


11.8 




fmax (kms -1 ) 


241 




r s (kpc) 


6.6 




<^red 


1.1 



a as fit parameters, and take = 3 kpc and ri j = 
0.3 kpc in rough accordance with observation and other 
mass models (e.g. [29 ). 

For the CDM(IT) model, we use M M , M ljb , v max 
and r s as fit parameters, and take ri jd - and ?r ; t> as fixed 
parameters with the MDM values. With data out to 
2i? (where Rq = 8.5 kpc) from [TT], [S] and [T3], we 
are able to generate good rotational curve fits for both 
models (see Figure [5] and Table [5]). 



600 - 



400 - 



200 - 



10 



15 
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Fig. 6 Predicted MDM (solid) and CDM(IT) (dashed) es- 
cape speed curves out to 2i£ for stars in the Galactic 
plane, using the parameters given in Tab. 2. The RAVE 
value 
Rq = 



u ESC = 544_ 46 kms 1 for the local escape speed (at 
8.5 kpc) is included as a reference point. 



Although the fits are simplistic, they are realistic 
enough to serve as adequate models of the Galaxy. The 
total visible (stellar) mass for the MDM fit is 6.1 x 
10 10 M Q with a bulge contribution of 1.3 x 10 10 M ; 
these values are comparable to estimates by [20], [17], 
[3"6] and [55]. The dark-to- visible mass ratio is /3 = 12, 
which approaches the cosmological ratio p%/p\ w 10 
in support of the hypothesis by BPR. The local sur- 
face mass density (at Rq — 8.5 kpc) of the disc is 
50 AT pc~ 2 , in good agreement with estimates by [2"1] . 
US], [7] and PI]. 

The total visible mass for the CDM(IT) fit is 6.6 x 
10 10 M Q with a bulge contribution of 1.2 x 10 10 M , 
while the local surface mass density of the disc is 56 M 
pc~ 2 . Again, these values match up reasonably well 
with other estimates. The local dark matter density is 
0.35 GeVcm -3 , which is towards the low end of bounds 
given by [T3] and [25]. 

Using the parameters given in Table [2J we obtain 
MDM and CDM(IT) escape speed curves for the Galaxy 
as in Section[3l but with a bulge term added to the grav- 
itational potential (see Figure©. For the MDM model, 
'/'bulge (f) is obtained in identical fashion to Equation 
([IB]) by numerically integrating Equations (fTo| and p7]l 
over the bulge; for the CDM(IT) model, it is given by 
Equation [fT5]), 

Again, escape speeds for stars in the Galactic plane 
are higher in the MDM model, with an offset that grows 
from « 20 kms" 1 to 50 kms -1 but remains nearly 
constant beyond r — Rq . The relative offset at r = Rq 
is Z\w/-ucdm(it) ~ 0.09. Both models predict values for 
the local escape speed that fall just short of the lower 
bound on the RAVE value v csc = 544+4g kms -1 ; in the 
CDM(IT) case, the fairly large discrepancy is due to 
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the truncation of the dark halo at r = 2R & . Indeed, by 
extending the truncation radius of the halo to 45 kpc, 
we find an escape speed for the CDM(IT) model that 
matches the observed RAVE value. Regardless, the fact 
that the MDM model is able to predict a local value of 
v CBC significantly greater than y/2v c i T corroborates the 
RAVE measurement, and lends credence to its own vi- 
ability as a dark matter model. 

5 Discussion and conclusions 

The addition of a non-luminous component of matter 
distributed spherically about the centre of a galaxy and 
extending well beyond the luminous disc component 
(i.e. a dark matter halo) explained the flat rotational 
curves of disc galaxies in the framework of Newtonian 
gravity, and led to the CDM paradigm in cosmology. In 
these models, a significant fraction of the gravitating 
matter in the Universe is in a dark collisionless form, 
with a density parameter ratio between dark and bary- 
onic matter of fid/ fib ~ 5. The CDM model has been 
extremely successful in explaining early Universe ther- 
modynamics, through to fluctuations in the CMB and 
the formation of structures on cosmological scales. Any 
viable alternative model must not only explain galactic 
rotational curves, but also match the successes of CDM 
in large-scale cosmology. 

A much-discussed alternative for explaining galaxy 
rotational curves is the modified Newtonian dynamics 
(MOND) model of [2&\, which entertains the possibility 
of an acceleration due to gravity that departs from the 
Newtonian form in regimes where the dynamical accel- 
erations are small. An in-depth review of the theory is 
provided by [33) . First proposed as a phenomenological 
model, the essence of MOND has subsequently been 
incorporated in a theory of modified gravity (MOG) 
by [3D], which is based on a covariant generalisation of 
Einstein's theory with the addition of auxiliary grav- 
itational fields to the metric. MOG has been success- 
fully used to model a large sample of galactic rotational 
curves [9] and galaxy clusters [10] ■ However, these mod- 
els that focus on the gravitational force appear not to 
be as well-placed to address questions on early Universe 
thermodynamics and large-scale cosmology. 

On the other hand, the MDM model of BPR is based 
on a fundamental theory of all particles and forces; as 
with the CDM model, it appears to have the poten- 
tial to not only explain galaxy rotational curves, but 
also to address problems relating to the early history 
of the Universe. In proposing their model, BPR point 
out that to date there has been no observational evi- 
dence for the dark matter postulated in CDM models. 
They note that one of the major unanswered questions 



in CDM cosmology, namely the reason for the close- 
ness of the observed ratio fid/ 'fib to unity, may find an 
easier solution in MDM models — where the dark (mir- 
ror) and baryonic (ordinary) sectors are expected to be 
closely linked. They also note that the density profiles 
predicted by n-body simulations for the CDM halos in 
galaxies show a central cusp, for which there is no obser- 
vational evidence from rotational curves of low surface 
brightness galaxies. While acknowledging that this may 
be a problem with the numerical simulations of galaxy 
formation in the CDM scenario, they point out that 
dark matter is collisional in the MDM model, and such 
cusps are not expected. 

Rotational curves, being based on the circular speed 
(usually of gas and sometimes augmented by that of 
stars) in galaxies, cannot by themselves provide strong 
constraints on dark matter models or gravity models. 
Indeed, CDM, MDM and MOG models have all been 
successfully used to model rotational curves of both low 
and high surface brightness galaxies. In this paper, we 
have sought to highlight some of the differences that 
may be expected in the dynamics of stars in galactic 
discs, which may potentially be used to discriminate 
between the CDM and MDM models. In particular, we 
have focused on stellar escape speeds from the galactic 
disc. Following BPR, our basic premise is that due to 
the collisional and dissipative nature of MDM, matter 
in the two sectors will have similar density distribu- 
tions up to a scaling factor. A consequence is that the 
MDM model has a very definite prediction on escape 
speeds that is readily falsifiable. Notwithstanding the 
uncertainties inherent in modelling the structure and 
extent of the CDM halos, we predict systematic and 
measurable differences in escape speed between the two 
models. 

We have taken the galaxies NGC 2403, DDO 170 
and DDO 154 (whose rotational curves have been mod- 
elled by BPR using MDM) along with a couple of syn- 
thetic galaxies as representative of spiral galaxies with a 
range of sizes, and used them to calculate and compare 
stellar escape speeds between the MDM and CDM (IT) 
models. We find that the escape speeds from the galac- 
tic disc are typically higher in the MDM model, with the 
percentage increase in escape speed relative to what is 
expected from an equivalent CDM(IT) model approach- 
ing Av/vcbm > 20 percent in the smaller low sur- 
face brightness galaxies. A similar procedure has been 
adopted for the Milky Way, where in addition to the 
MDM stellar disc, we have also included an MDM bulge 
component. 

In the case of the Milky Way, recent estimates from 
the RAVE survey yield an escape speed from the so- 
lar neighbourhood of v csc = 5441^ kms -1 , and an 
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escape to circular speed ratio v csc /v c i r ~ 2.5. If one 
assumes normal gravity and axially symmetric density 
structures, the circular speed gives a measure of the 
amount of matter within the solar circle, while the es- 
cape speed gives a measure of the amount of matter 
outside the solar circle. A measurement of v esc /v c i T that 
is significantly larger than y/2 could then be taken as 
independent evidence that there is a large proportion of 
mass outside the solar circle. However, we have shown 
that our MDM model of the Milky Way yields a simi- 
lar value Vosc/^cir ~ 2.2, and that given the errors, this 
particular test cannot be used to establish the presence 
of a dark halo or to distinguish between the two models 
at the present time. 

This situation is expected to change when the Gaia- 
ESO public spectroscopic survey using FLAMES be- 
comes available in 2016 [2T]; the survey will yield accu- 
rate position and velocity information of a large number 
of individual stars in the thin and thick discs, the bulge 
and the halo. Detailed studies of the motions of stars in 
the solar neighbourhood (as has been done with RAVE) 
will allow a more precise measurement of the local es- 
cape speed. With the new set of data, the study could 
be extended to establish the radial dependence of the 
escape velocity, which will provide even stronger con- 
straints on the different dark matter models. We may 
expect substantial differences under the MDM model in 
the kinematics of stars in the galactic halo, which may 
also be measurable in future surveys. 

The measurement of escape speeds from discs in 
external galaxies may also soon become possible. This 
will require the accumulation of position and velocity 
data for a large number of individual stars using multi- 
object optical and infrared spectroscopy with large op- 
tical telescopes, along the lines of the recent spectro- 
scopic study of the nearby dwarf irregular galaxy WLM 
using the VLT and Keck II [35]. Statistical studies of 
the resulting stellar velocity distributions of a number 
of such galaxies viewed at different inclinations could 
then be used to estimate escape speeds and to test the 
MDM paradigm. 
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